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ABSTRACT 

We present NGC 6221 as a case study of "X-ray loud composite galaxies," which appear similar to 
star bursts at optical wavelengths and resemble traditional active galactic nuclei in X-rays. The net 
optical spectrum of NGC 6221 is properly characterized as a starburst galaxy, but in X-rays, NGC 
6221 is similar to Seyfert 1 galaxies, exhibiting a power-law continuum spectrum, a broad Fe Ka line, 
and continuum variability on timescales of days and years. High-resolution images reveal that the 
detected active nucleus is relatively weak, not only at optical, but also at near-infrared wavelengths. 
An obscuring starburst, in which the interstellar gas and dust associated with the starburst conceal the 
active nucleus, accounts for these peculiar features. We demonstrate quantitatively that obscuration by 
column density Nh = 10^^ cm~^ combined with relatively weak intrinsic nuclear activity can produce 
an optical spectrum characteristic of the surrounding starburst alone. While optical surveys would not 
identify the active nuclei that make these galaxies significant X-ray sources, such galaxies may in fact 
be important contributors to the X-ray background. 

Subject headings: Galaxies: individual (NGC 6221) — galaxies: Seyfert — X-rays: galaxies 



1. INTRODUCTION 

Since the discovery of the cosmic X-ray background 
(XRB) (Giacconi et al. 1962), great progress has been 
made in understanding its origin. Surveys with ROSAT 
successfully resolved 70-80% of the soft X-ray background 
at energies 0.5 to 2 keV into discrete sources (Hasinger 
et al. 1998; Schmidt et al. 1998). Current work with 
the Chandra X-ray Observatory demonstrates that dis- 
crete sources account for similar fractions of the harder 
X-ray background, at energies up to 10 keV (Giacconi et 
al. 2001). The problem is not completely solved, however. 
Although much of the XRB is made up of what we might 
classify as the nuclei of otherwise normal bright galaxies or 
typical active galactic nuclei (AGNs), a significant fraction 
of the discrete sources are not easily identified. If they are 
AGN, they appear to have unusual properties. Mushotzky 
et al. (2000) find a population of optically faint sources, 
which could either be the active nuclei of dust-enshrouded 
galaxies or the first quasars at very high redshifts. The 
obscuration of some AGNs is so strong that they are not 
detected in X-rays at all. If the SCUBA submillimeter 
sources, for example, are not driven entirely by star for- 
mation, they may contain AGNs with Compton-thick tori 
and little circumnuclear X-ray scattering (Hornschemeier 
et al. 2000). 

Among nearby X-ray emitters are galaxies that have 
X-ray characteristics of AGN, but whose optical spectra 
are typical of H ii regions or starbursts. These objects 
may comprise some of the "unusual" discrete sources of 
the XRB, seeming innocuous in larger optical surveys, 
but generating significant hard X-ray emission. Moran, 
Halpern, & Helfand (1996) dub them "composite galaxies" 



because they share some characteristics of both starbursts 
and AGNs, albeit at different energies. We prefer "X-ray 
loud composite galaxies," however, to distinguish these ob- 
jects from the larger class of AGN/starburst composite 
galaxies that exhibit their dual nature consistently at all 
wavelengths (e.g., Levenson, Weaver, & Heckman 2001a; 
Cid Fernandes et al. 2001). While physically, members of 
both classes contain AGN and starburst components, their 
observational characteristics are distinct. 

The X-ray loud composite galaxies comprise a small 
fraction of soft X-ray selected sources. Moran et al. (1996) 
examined optical spectra of the Boiler et al. (1992) catalog 
of IRAS sources detected in the ROSAT AlLSky Survey 
and found 7 of this type among the 210 galaxies they clas- 
sified, and they suggest that obscuration, including of the 
narrow-line region, may play a role in accounting for their 
unusual properties. The optical spectra have the charac- 
teristics of starburst and H ii region systems, rather than 
AGNs, based on diagnostic diagrams of their emission line 
ratios (e.g., Veilleux & Osterbrock 1987). Only on closer 
examination do the optical spectra hint at additional ac- 
tivity, such as broadened [O ill] lines or a weak broad com- 
ponent underneath the Balmer lines. At X-ray energies, 
however, the X-ray loud composites appear to be typical 
AGNs, sharing characteristics of luminosity and variabil- 
ity. 

Here we present NGC 6221 as a case study of the X- 
ray loud composite galaxies. Advantageously, this galaxy 
has been observed over a wide energy range, and we uti- 
lize ground-based optical spectroscopy, high-resolution op- 
tical and near-infrared images from the Hubble Space Tele- 
scope (HST), and X-ray imaging and spectroscopy from 
ASCA and ROSAT to physically account for the emis- 
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sion sources. We adopt a distance of 19.8 Mpc to NGC 
6221, with Hq = 75kms~^ Mpc^^ for recession velocity 
cz = 1482kms-\ so 1" = 96 pc. Although NGC 6221 has 
been classified as a Seyfert 2, its optical spectrum is more 
characteristic of a starburst galaxy. The X-ray emission of 
NGC 6221 certainly qualifies it as an active galaxy, and, 
as we argue below, it is more like a type 1, in fact. We 
propose the "obscuring starburst" model to account for 
these properties and demonstrate in general that a star- 
burst can easily hide an AGN. While previous studies of 
the starburst- AGN connection have faced the initial task 
of finding starbursts in active galaxies (e.g., Heckman et 
al. 1997; Gonzalez Delgado et al. 1998; Gonzalez Delgado, 
Heckman, & Leitherer 2001), the X-ray loud composite 
galaxies represent the other extreme composition, in which 
the dominant starburst conceals the weaker AGN at opti- 
cal wavelengths. 

2. OPTICAL SPECTROSCOPY 

2.1. Observations 

A long-slit optical spectrum of NGC 6221 was obtained 
on the night of 1997 September 24 with the CTIO 4-m 
telescope using the Cassegrain Spectrograph and Loral 3K 
CCD. The 1'.'5- (144 pc-) wide slit was oriented along the 
paralactic angle at position angle 68°. The spectral range 
covered was 3650-7300A, at a resolution of approximately 
4A. The seeing was 1'.'2. After the long-slit spectrum was 
reduced and flux-calibrated, one-dimensional spectra were 
extracted. We binned together 3 pixels along the slit, 
which corresponds to 1'.'5 or 144 pc, except the outermost 
spectra, centered 5'.'25 from the nucleus, in which 6-pixel 
(3" = 288 pc) regions were combined. 

2.2. The AGN Narrow-Line Region 

The net optical spectrum of NGC 6221 is typical of a 
reddened starburst (Phillips 1979; Morris & Ward 1988; 
Storchi-Bcrgmann, Kinney, & Challis 1995). The only 
sign of non-stellar activity at these wavelengths is an 
[O III]A5007 component broader than and blue-shifted 
with respect to H/3. This feature, as well as the early detec- 
tion of NGC 6221 as an X-ray source (Marshall et al. 1979) , 
motivated Veron, Veron, & Zuiderwijk (1981) to propose 
a composite Seyfert 2/star burst scenario, in which a faint 
AGN narrow-line region (NLR) is superposed on a dom- 
inant H II region powered by massive stars. Our spectra 
confirm the presence of the NLR-Iike [O III] component, as 
illustrated in Figure 1, where we show the [O III] and H/3 
profiles for the nucleus and six off-nuclear positions along 
the slit. In the nuclear (1'.'5 x 1'.'5) spectrum, the broad 
and blue component dominates the [O III] profile, extend- 
ing to more than 1000 kms~^, while in H/3, the narrow 
component dominates. 

The flux ratio, R = F[oiii]/Fhi3, increases from 0.3-0.5 
(typical of starburst nuclei) within 100 kms^^ of the H/3 
peak to ii > 2 (typical of Seyferts) for v > 500 kms~^ 
toward the blue, in agreement with the results of Veron et 
al. (1981) and Pence & Blackman (1984). We agree with 
previous studies that the broadened feature likely origi- 
nates in an AGN NLR. The observed high R in the blue 
wing of the lines can in principle be achieved in metal poor 
starbursts, but this would be inconsistent with the approx- 



imately solar mctallicity inferred for NGC 6221 (Storchi- 
Bergmann, Calzetti, & Kinney 1994; see also Durret & 
Bergeron 1988). The width of this high excitation compo- 
nent as well as its flux ratio is more typical of AGNs than 
of starburst galaxies. 




V [km s"'] 

Fig. 1. — H/3 (solid line) and [O III] {dotted line) line profiles 
across the slit. The five central extractions are I'.'b wide, while each 
of the extreme extractions averages the emission from 3'.' 75 to 6'.' 75 
from the nucleus. Profiles are normalized from to 1 and offset 
by constant values for ease of comparison. The actual [O III] /H/3 
integrated flux ratio ranges from 0.5 to 0.9. The angular scale for 
NGC 6221 is 96pcarcsec~^, so the profiles shown here correspond 
to distances up to 500 pc from the nucleus. 

We quantify the characteristics of the two components, 
modelling the line profiles with two Gaussians. We per- 
formed a global fit of the H/3 and Ha -|- [N II] profiles using 
common kinematic parameters (centroid and FWHM) for 
all lines. The [S II] lines were not used in this fit to deter- 
mine the kinematic parameters but were later deblended 
with the two resulting Gaussians. In all extractions, these 
combined fits yield a broad component of FWHM ^ 500- 
600 kms^^ (after correcting for instrumental broaden- 
ing) blue-shifted by 150-250 kms^^ with respect to a 
narrower unresolved component (whose intrinsic FWHM 
~ 100 kms"^; Veron et al. 1981). The [O III] line was 
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modelled separately, with one component fixed to have the 
same centroid and FWHM of the narrow Gaussian identi- 
fied in the H/3 + Ha + [N II] fits and the other component 
unconstrained. This yields a more satisfactory fit to the 
[O III] profile because, even though all lines have the same 
general structure, the broad component is much stronger 
in [O III]. Fitting the [O III] yields a broad component 
with FWHM - 600 kms"! shifted by approximately 230 
kms^^ with respect to the narrow unresolved component, 
which better represents the kinematic properties of the 
NLR. 

With this evidence of line width and flux ratios, we refer 

to the broad component as the NLR or AGN component 
and the narrow one as the starburst or H II component, 
corresponding to gas photoionized by the AGN and mas- 
sive stars, respectively. Table 1 lists the fluxes of the main 
emission lines for flve positions along the slit. The per- 
centage of the flux contained in the NLR component is 
listed in parentheses except in the case of [O I], to which 
we have not applied the profile decomposition. 
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Emission Lines 
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l'.'5NE 
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l'.'5SW 
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H/3 


100 (23) 


100 (29) 


100 (37) 


100 (35) 


100 (27) 
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Ha 
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700 (30) 
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690 (36) 


660 (21) 


[N II]A6584 


420 (30) 


380 (32) 


420 (40) 


360 (36) 


400 (17) 


[S II]A6717 


79 (25) 


63 (33) 


73 (45) 


83 (45) 


95 (26) 


[S II]A6731 


76 (29) 


64 (31) 


82 (41) 


82 (41) 
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Fif(9 (erg cm-^ li- 


') 1.5 X 10-" 


1.1 X 10-" 


2.4 X 10-'" 


1.2 X 10-" 


4.0 X 10-" 


ne (cm-^) 


500 


670 


710 


580 


220 


Av (mag) 


2.8 


2.8 


3.0 


2.8 


2.6 



spatially resolved. The [O III] asymmetry is easily detected 
in the extractions centered 3" (corresponding to 290 pc) 
away from the nucleus, and we marginally detect it in the 
outer (and noisier) spectra. This projected size of 300- 
500 pc is typical of observed NLR sizes in Seyferts (e.g., 
Schmitt & Kinney 1996). While the blue [O III] wing grad- 
ually fades farther from the nucleus, its profile increasingly 
resembles that of H/3, and the NLR share of the total line 
flux decreases (Table 1). 

The fact that the NLR feature is seen on both sides of 
the nucleus gives an indirect clue of the geometry of the 
AGN-photoionized gas in NGC 6221. This symmetry im- 
plies that the NLR is not projected as a one-sided cone. 
More importantly, the broad [O III] component is also ex- 
tended along several other position angles (H. A. Fraquelli 
& T. Storchi-Bergmann, in preparation), so the NLR is 
approximately symmetric in NGC 6221. If the intrinsic 
shape of the NLR is roughly conical, then these results in- 
dicate that we see it from within the opening angle of the 
cone. Projected onto the plane of the sky, the NLR thus 
appears nearly circular. This geometry implies that if we 
could turn off the starburst, NGC 6221 would be classified 
as a Seyfert 1. 

2.3. No Sign of a Broad-Line Region 

There is no indication at all of broad (FWHM > 1000 
kms~^) permitted emission lines in NGC 6221, character- 
istic of Seyfert Is, even at longer wavelengths (Moorwood 
& Oliva 1988). (We stress that the "broad" component we 
refer to above is not associated with an AGN Broad-Line 
Region (BLR), since (1) it is present in [O III], and (2) it is 
too narrow for BLR standards.) Given that the only mea- 
surable optical signature of an AGN in NGC 6221 is the 
relatively subtle broadened [O III] feature, we determine 
the upper limit on BLR characteristics in the optical spec- 
trum. By adding a BLR-like (FWHM between 2500 and 
7000 kms"-'^) component to the observed H/3 profile, we 
estimate that any such BLR component has a maximum 
equivalent width EWf of 10 A; in comparison, the mea- 
sured equivalent width of H/3 is 21.5 A. In terms of flux, 
this limit corresponds to F^^^ < 1.1 x 10~^^ergcm~^ s~^. 



Note. — Line fluxes relative to H/3 for five 1'.'5 X 1'.'5 extractions along the slit described in the 
text. Numbers in parentheses are the percentage of the line flux that is contained in the "NLR 
component" identified in the profile fitting. 

All the total nuclear emission line ratios are consistent 
with a starburst classification on conventional diagnostic 
diagrams: [O III]/H/3 = 0.89; [N II]A6584/Ha = 0.56; 
[O I]A6300/Hq; = 0.02; and [S II]AA6717 -h 6731/Ha = 
0.21. We note, however, that [N II] /Ha lies at the bound- 
ary between AGN and H Il-region-like objects when plot- 
ted against [O III]/H/? (Veilleux & Osterbrock 1987). 
While the [O III] /H/3 ratio for the broad and narrow com- 
ponents yields NLR and H Il-region-like ratios (2.1 and 
0.2 respectively), the [N II]/Ha and [S II]/Ha ratios are 
only slightly larger for the broad component than for the 
narrow one. Because [O III] is the brightest line in NLRs, 
it is not surprising that our simple profile decomposition 
does not disentangle NLR and H Il-like line ratios for [N II] 
and [S II] as effectively as for [O III] . Our modelling of the 
starburst and NLR components in §5 uses only the results 
for the H/3 and [O III] lines, so we have not attempted a 
more detailed profile decomposition. Table 1 also lists the 
electron density inferred from the total [S II] line emission, 
calciflated using the IRAF'' task temden. The density de- 
creases outwards, from rig ^ 700 cm""^ in the nucleus to 
100-300 cm~^ at about 500 pc, in common with other 
starburst systems (Lehnert & Heckman 1996). 

Besides revealing the NLR component of the AGN, Fig- 
ure 1 shows that the region producing this component is 

IRAF is distributed by the National Optical Astronomy Observatories, operated by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National Science Foundation. 



which, correcting for Ay = 3 mag (see below), translates 
into an upper limit Lf/fj^ < 1.3 x 10^*^ erg s"-"^ on the cor- 
responding intrinsic luminosity. 

A Seyfert galaxy without detectable broad optical emis- 
sion lines would be classified as a type 2, as Veron et 
al. (1981) originally proposed for NGC 6221, but several 
pieces of evidence indicate that its orientation is charac- 
teristic of Seyfert Is. First, the extended NLR emission 
noted above suggests that we view the nuclear region face- 
on. Second, all the X-ray properties, including luminos- 
ity, variability, and emission line attributes, are typical 
of Seyfert Is. We return to this issue in §5, developing a 
physical model of the geometry of NGC 6221 that explains 
the absence of broad optical lines despite the orientation. 

2.4. Extinction 

A distinctive feature of NGC 6221 is its large extinc- 
tion. The values of Ay listed in Table 1 were computed 
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Table 2. Summary of Emission Line Properties 
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Note. — NLR and siarburst (SB) components of 1 he main emission lines for the nuclear and integrated spectrum. De-reddened 
quantities are corrected for Ay — 3 mag. Fluxes are in units of ergcm-^ s-i and luminosities in ergs-i. 



from the measured Ha/H/3 ratio using the extinction law 
of Cardclh, Clayton & Mathis (1989, with Ry = 3.1) and 
assuming a case B intrinsic ratio of 2.86 (Osterbrock 1989) 
appropriate to starburst nebulae. In the nuclear extrac- 
tion, we measure Ay = 3.0 mag, compatible with the esti- 
mate of Storchi-Bcrgmann ct al. (1994), and wc agree with 
their conclusion that the correction for underlying absorp- 
tion in the Balmer lines has a negligible effect on the es- 
timated extinction, so we have not considered it in our 
analysis. The Ha/H/? ratio for the NLR component is 7.9, 
somewhat larger than the 7.2 obtained for the starburst 
component. This suggests that that the NLR suffers an 
additional 0.3 mag extinction with respect to the starburst 
lines. Since this is a small difference, we adopt a common 
visual extinction of 3 mag for the entire emission line re- 
gion in NGC 6221. Dust within the Milky Way accounts 
for Ay = 0.5 mag (Schlcgcl, Finkbcincr, & Davis 1998), so 
the intrinsic extinction in NGC 6221 Ay = 2.5. We sum- 
marize the observed fluxes, de-reddened fluxes, and lumi- 
nosities of the NLR and starburst components of several 
important lines in Table 2 for both the nuclear extraction 
and the total spectrum integrated from — 6'.'75 to -|-6'.'75 
along the slit. 

2.5. Properties of the Starburst in NGC 6221 

We characterize the starburst in NGC 6221 using the 
Starburst 99 models of Leitherer et al. (1999). For a 

constant star formation rate model with solar metallic- 
ity and Salpeter IMF between 1 and 100 M0, the star 
formation rate is related to the Ha luminosity as SFR = 
-f'Ha/(2.98 X lO'^^ ergs~^)MQ yr"^ Applying this relation 
to the starburst component in the nuclear spectrum (Ta- 
ble 2) gives SFR = O.17Af0 yr"! or 8 x IQ-^M© yr"! pc-^. 
Having such a high central SFR, NGC 6221 certainly quali- 
fies as a true starburst, rather than a normal galaxy, com- 
pared with the samples of Kennicutt (1998a), for exam- 
ple. The strongest star formation is concentrated within 
the central region, extending to a radius of approximately 
400-500 pc, and it is not confused with ongoing star for- 
mation in the galaxy's disk. 

Within the limited circumnuclear region of NGC 6221, 
however, the star formation is widespread, as the spec- 
tral extractions and HST images (§3) demonstrate. This 
can also be seen in Figure 43 of Cid Fernandes, Storchi- 
Bergmann, & Schmitt (1998), where the equivalent widths 
of absorption lines and continuum colors are plotted as a 
function of distance from the nucleus. The bluer colors and 
highly diluted absorption lines in the extended central re- 
gion are typical of other starburst and Seyfert 2/starburst 
composite galaxies observed by those authors and others 



(Gonzalez Delgado et al. 2001). The widespread star for- 
mation is further confirmed by Storchi-Bcrgmann et al. 
(1995), who measure the Ha flux in a 10" x 20" aperture 
to be 5.4 times larger than the Ha flux we measure in our 
1"5 X 1"5 nuclear window. 

The far-infrared (FIR) luminosity of NGC 6221 Lfir = 
2.7 X 1O^"L0. Using the bolometric correction of Meurer, 
Heckman, & Calzetti (1999), this translates into Li,oi = 
3.7 X W^^Lq. a constant SFR model predicts Lboi ~ 
1.1 X per A^Qyr"^ of star formation, so the im- 

plied total SFR = 3.3 M0yr~^, 20 times larger than the 
rate inferred for the nuclear (144 x 144 pc^) region alone. 

3. HIGH-RESOLUTION HST IMAGES 

HST acquired optical and near-infrared (NIR) images 
of NGC 6221, which we obtained from the HST archive 
at the Space Telescope Science Institute. The Wide-Field 
Planetary Camera 2 observation provides spatial resolu- 
tion of about 0"l(w 10 pc) in the 0"046 pixels of the Plane- 
tary Camera (PC), which covered the nuclear region. The 
F606W filter was in place, for broad coverage at optical 
wavelengths; the total system throughput is centered at 
wavelength = 6030A, with width 5X = 1500A. We use 
the standard pipeline processing of this single 500 s ex- 
posure from 1995 April. The Near Infrared Camera and 
Multi Object Spectrometer (NICMOS) camera 2 observed 
NGC 6221 for 320 s in 1998 May with the F160W filter 
(Ac = 1.60/im, (5A = 0.40/xm). The camera's plate scale is 
0'.'075 pixel- \ for 0'.'2 resolution. 

We aligned the two images in order to identify the galaxy 
nucleus and other sources common to both, although sub- 
sequent photometry was performed on the original cali- 
brated data. First we used the nominal astrometry to 
perform a two-dimensional polynomial transformation of 
the NICMOS image to the PC orientation and scale. Sec- 
ond, we measured the centroids of five point sources com- 
mon to both frames and located outside the nuclear region. 
We then linearly shifted the transformed NICMOS frame 
along the rows and columns by the average measured off- 
sets. 

We display the results in Figure 2. We identify the 
brightest central NIR source as the AGN, with colors and 
luminosity similar to other active galaxies, as we discuss 
below. The source is evident in the optical image, but it is 
not the brightest one. In the optical band, the nucleus is 
relatively weak compared with the diffuse emission and the 
other bright sources, which we identify as stellar clusters. 

We measure the brightness of the nuclear source alone 
using point spread function (PSF) subtraction. Because 
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Fig. 2. — Co-aligncd images of the eenter of NGC 6221 in the optical (left) and near-infrared {right). The brightest central source detected 
at 1.6jUm is the AGN, although it is relatively weak at optical wavelengths. 



the PSF varies significantly across the field of the PC, 
we obtained a calibration observation of an unsaturated 
point source located at a similar position in the PC. We 
scaled this empirical PSF to match the flux of the AGN 
within the central 3x3 pixel background-subtracted re- 
gion, computing the total fiux from the integration of the 
entire scaled PSF. In the F606W filter, which is compara- 
ble to a F filter, the net result is meoeo = 18.1 mag on the 
Vega scale, or flux density /eoeo = 0.19 mJy. Because the 
NICMOS PSF does not vary significantly across the field, 
we use a point source within the field of view to define 
the PSF. The PSF fitting yields mis = 13.6 mag on the 
Vega scale in the F160W filter, which is comparable to the 
standard H band, or /i.g = 3.8 mJy. At 1.6/im, the PSF 
subtraction reveals an incomplete ring of diffuse emission 
whose radius is approximately 0'.'3 and hints at a source 
about Cf.'l away from the nucleus, which is clearly resolved 
in the optical image. 

The central point source is very red compared with 
quasars, and therefore must suffer substantial extinction 
if it is the active nucleus. Specifically, we compare with 
the V — H colors of quasars that Elvis et al. (1994) ob- 
served. Although this sample includes radio-loud as well 
as radio-quiet objects, it is appropriate for comparison be- 
cause the relative contribution of the host galaxies is small, 
so it yields reasonable measurements of the colors of the 
nuclei alone. The nucleus of NGC 6221 equals the median 
V — H = 2.1 mag of the sample, assuming standard inter- 
stellar extinction {Ah /Ay = 0.19; CardeUi et al. 1989) for 
Ay = 3.0, the same value deduced from the steep Balmer 
decrement measured in the optical spectrum (§2.4). 

The high-resolution images clearly reveal that the nu- 
cleus of NGC 6221 is weak with respect to the surrounding 
starburst, especially at optical wavelengths. Through the 
F606W filter, the nucleus provides only 6% of the observed 
flux within the 1'.'5 x 1'.'5 aperture of the central region in 
which the central optical spectrum was extracted (§2), as 
Figure 3 graphically illustrates. Even at 1.6/im, where the 
AGN is relatively bright, it accounts for only 32% of the 
flux within a 1'.'5 x 1'.'5 region. 
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Fig. 3.— Central region of NGC 6221 observed with the HST 
PC. A circle of O^'ll radius, corresponding to the width of the point 
spread function, is drawn around the active nucleus. The overlaid 
square is 1'.'5 on a side, which corresponds to the aperture of the 
central optical spectrum. 

We calculate H/? flux of the AGN broad- line region from 
these measurements of the point source. For / oc z^", the 
spectral slope between the H and V bands a = —3.1, so 
/4861 = 0.096 mJy. The equivalent width of H/3 in active 
galaxies is approximately lOOA (e.g., Boroson & Green 
1992; Binette, Fosbury, & Parker 1993), nearly aU of which 
is due to the BLR in type 1 sources. Thus, we estimate the 
observed (and extincted) F^^^ = 1.2 x lO"^** erg cm^^ 
from the combined HST data. This value is at the de- 
tection limit of the optical spectrum and higher than the 
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NGC 6215 



NGC 6215A 



NGC6221 



5' 




8-GOl 



Fig. 4. — NGC 6221 field observed with ASCA GIS2. North is at the top, and east is toward the left. The minimum contour is 2cr above 
the background, and levels axe scaled logarithmically by factors of 2. The background is not subtracted from the greyscale image, to illustrate 
the circular field of view. 



estimates of the intrinsic AGN strength based on X-ray lu- 
minosity we compute below (§5). Note that although the 
extinction in the lines and continuum may be somewhat 
different (e.g., Witt, Thronson, & Capuano f992), in this 
rough calculation, wc have reasonably assumed that they 
are equal. 

While the evidence for identifying this red central source 
as the AGN is indirect, it offers the most plausible account 
of all the data. First, the observed color of this source is 
typical of active nuclei, given the reddening we measure in 
the spectrum. Second, the V, H, and X-ray luminosities 
simultaneously match the spectral energy distribution of 
a typical quasar (§5). Third, the expected line emission 
of the point source is consistent with the upper limit on 
detection of the BLR. Finally, none of the other detected 
sources share these characteristics, so they are unlikely to 



be the active nucleus. We could fail to detect the AGN 
entirely, even at 1.6/um. Within the bright central region 
of NGC 6221, we would not detect a source fainter than 
mi. 6 = 18 on the Vega scale. In this case, the intrin- 
sic AGN would be extremely weak. While a faint, buried 
AGN is certainly possible, it would leave the problem of 
explaining the nature of the central source, with its excep- 
tionally red color compared with the surrounding stellar 
clusters and diffuse emission. Near-infrared spectroscopy 
of the red central knot would clearly test its identifica- 
tion as the AGN. High spatial resolution is feasible, and 
the contrast between this source and its surroundings is 
greater at NIR than optical wavelengths. The NIR char- 
acteristics of AGN, including high-ionization lines and the 
BLR component of hydrogen recombination lines, could 
be observed. 
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4. X-RAY IMAGING, SPECTROSCOPY, AND VARIABILITY 

4.1. Data Reduction 

NGC 6221 was observed with the ROSAT HRI for 12 

ks in 1995 September and with ASCA in 1997 September, 
resulting in exposure times of 35 ks and 37 ks in the SIS 
and GIS detectors, respectively. The ASCA data were pro- 
cessed to exclude times of high background rates and to be 
appropriate for the SIS 1-CCD observation mode. The ac- 
cepted data were selected to be outside of and greater than 
16 s after passages through the South Atlantic Anomaly, 
at elevations above the Earth's limb greater than 5°, at 
times greater than 16 s after a satellite day/night tran- 
sition, and with geomagnetic cutoff rigidity greater than 
6 Gev c^^ . We grouped the spectra into bins having a min- 
imum of 30 counts, so statistics are applicable. 

4.2. NGC 6221 

The X-ray data are described in detail in Levenson, 
Weaver, & Heckman (2001b), and we summarize these 
results below. The ASCA image of NGC 6221 obtained 
with the GIS2 detector (Figure 4) illustrates the X-ray- 
emitting environment of this galaxy. NGC 6221 is in a 
group. Extended diffuse emission around NGC 6221 and 
ESQ 138-GOl implies that these galaxies are dynamically 
related, and NGC 6215A may also be included in this 
group. Observations at 21-cm and a similar recession ve- 
locity indicate that NGC 6215 may be interacting with 
NGC 6221 (Pence & Blackman 1984), possibly accounting 
for its peculiar morphology, although the X-ray emission 
from NGC 6215 is separated from the group's emission. 
The ASCA field of view also includes part of the nearby 
cluster 1RXSJ165259.4-594302, identified in the ROSAT 
All-Sky Survey (Voges et al. 1999). 

NGC 6221 is significantly extended in ROSAT YLm ob- 
servations, on a physical scale of at least 5 kpc. The re- 
solved component accounts for more than half the soft 
X-ray emission in this observation from 1995 September. 
By analogy with normal starburst galaxies lacking active 
nuclei, this extended component is likely due to heating of 
the interstellar medium by stellar winds and supernovae, 
possibly including a large-scale outflow where the starburst 
"superwind" has broken out of the galaxy's disk (Chevalier 
& Clegg 1985; Heckman, Lehnert, & Armus 1993). The 
characteristic X-ray spectrum of this diffuse component 
observed in pure starburst galaxies and active galaxies that 
are known to contain starbursts is soft thermal emission, 
which is due to gas that has been heated to temperatures 
of 0.2 to 1.0 keV. 

Spectrally, in the ASCA observations obtained two years 
later, NGC 622 1's AGN characteristics dominate, and the 
soft flux increased by a factor of 3. The total soft (0.5-2 
keV) X-ray luminosity Lo.5-2 = 1-4 x lO^'^ergs"^, and 
the hard (2-10 keV) X-ray luminosity -L2-10 = 6.6 x 
10^"'^ ergs^^, where both are corrected for Galactic absorp- 
tion and measured in the SISO detector. Levenson et al. 
(2001b) model the 0.5-10 keV spectrum with four com- 
ponents: 1) the AGN as a power law with photon index 
r = 1.9, absorbed by Nh = lO^^cm"^; 2) the scattered 
AGN spectrum, with 7% of the intensity of the intrinsic 
nucleus; 3) a broad (cr = 0.5 keV) emission line centered 
on = 6.6 keV; and 4) thermal emission with kT = 1.4 



keV. Because the X-ray absorption of the intrinsic AGN is 
small, even at softer energies (0.5-2 keV), 80% of the soft 
X-rays are due to the AGN viewed directly, and most of 
the remainder constitute the scattered component. Only 
5% of the soft emission included in the spectral model is 
thermal. 

The variation in X-ray intensity between the HRI and 
ASCA observations is due to changes in the AGN, a change 
of either its intrinsic luminosity or its obscuring column 
density. The measured flux difference requires cither an 
increase by a factor of about 5 in the intrinsic AGN lumi- 
nosity or a decrease in column density from Nh > 10^"^ at 
the time of the ROSAT observation to Nh = 10^^ cm'^. 
While the relative fraction of spatially extended X-ray 
emission, which we associate with thermal emission (Lev- 
enson et al. 2001b, a), diminished from more than half to 
only 5% of the soft flux, its total flux likely remained con- 
stant. 
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Fig. 5.— Light curves of NGC 6221 in GIS3 (top) and SISO {bot- 
tom) detectors. The count rates, extracted in 400-s bins for clarity, 
have been background-subtracted. The best-fitting constant rate 
for each case is overlaid. In neither case does the constant model fit 
the data acceptably. 

On shorter timescalcs, within the ASCA observation, 
we also measure significant variability. Background- 
subtracted lightcurves are shown in Figure 5 for the GIS3 
and SISO detectors, where 400-s bins clearly illustrate the 
significant variablity. The best-fitting constant count rate 
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Table 3. Other Sources in GIS Fields 



Object 


Distance 
(Mpc) 


GIS Rate 
(10-^ cts s-i) 


Nh 
(10^1 cm-2) 


r 


or kT'^ 






^^2-10 


^0.5-2 


^2-10 


NGC 6215 
NGC 6215A 
1RXSJ1652 


20.8 
38.6 


8.9 ± 0.97^ 
6.9±0.99f 
53 ±2.11! 


1.5 :t^} 
1.5f 


r = 
r = 

kT 


1 t;+0.9 

1 fi+°-^ 

^■"-0.6 
_ 9 O+0.9 


45/39 
42/46 
114/93 


rt 9 + 3.5 
^•^-1.2 
o 0+1.7 

09+3.8 
■'''-3.5 


-'- '-9.3 

7.4t|« 
4011^ 


-'-•-'--0.6 

0-^-2.6 





'^Power law photon index or bremsstrahlung temperature in keV. 
''0.5-2.0 keV model flux in units of 10~^^ ergcm"^ s~^. 
■^2.0-10.0 keV model flux in units of lO^^"* crgcm^^ s"^. 

'^O. 5-2.0 keV luminosity, corrected for Galactic absorption, in units of 10'^'' erg s^^. 
*'2. 0-10.0 keV luminosity, in units of lO^'^ erg s^'. 
'G1S3 count rate and fluxes. 
^G1S2 count rate and fluxes. 

Note. — Errors on model parameters are 90% confldence limits for two interesting parameters, except fluxes, where 
errors are 90% confidence for one parameter. Parameters that are constrained by hard limits (e.g.. Galactic column 
density) are marked with a colon. Fixed parameters are noted with f. 



is drawn, but this constant flux model fits the data un- 
acceptably in both cases. We extracted spectra over lim- 
ited time intervals to examine the (earlier) lower flux and 
(later) higher flux states separately. With the limited 
statistics of these shorter exposures, however, we could 
not identify significant spectral differences between the two 
states. The same physical model that fits the total spec- 
trum arc also appropriate for the time-constrained spectra, 
with measurable variation only in the intensity. Similarly, 
the hardness ratio, defined in terms of the hard and soft 
fluxes as (F2_io-^b.5-2)/(-F2-io+-Pb.5-2), shows no trend 
with the variation of total intensity. 

The X-ray flux varies by about 50% over a timcscalc of 
5 X 10^ s. We quantify the variation in terms of the "excess 
variance," u^ms) defined by Nandra et al. (1997a). Sam- 
pling in 256-s bins in order to compare with published re- 
sults, we find 0-2^, = 0.062 ±0.014 for the GIS3 lightcurve 
of NGC 6221. Thus, NGC 6221 is similar to broad-lined 
Seyfert 1 galaxies in the trend of increasing cr^ms with de- 
creasing hard X-ray luminosity; e.g., Turner et al. (1999) 
measure (Trms ~ 1 x 10"'' at ^2-10 ~ 10**^ ergs"-', rising 
to a^^, « 0.05 for ia-io « lO^^ergs-i. 

In addition to the unusual variability, its Fe line emis- 
sion makes NGC 6221 unlike most other Seyfert 2s. It 
contains a broad Fe Ka line centered at = 6.6 keV with 
EW- 400 cV and line width a = 0.5(+0.4, -0.02) keV 
with 90% confidence. The Fe line properties of this galaxy 
are more characteristic of Seyfert Is, where the reprocess- 
ing material is inferred to be close to the galaxy nucleus 
(Nandra et al. 1997b), as opposed to the narrow Ka lines 
at E — 6.4 kcV duo to neutral Fe farther from the nucleus, 
which is typical of Seyfert 2s, including Seyfert 2/star burst 
composites (Levenson et al. 2001b). 

4.3. Other Sources 

The group members NGC 6215 and NGC 6215A are in- 
cluded in the larger field of view of the GIS detectors, and 
we analyze their X-ray spectra here. NGC 6215 lies near 
the extreme edge of the field of view, so we extracted the 
spectrum from an elliptical region having semimajor axes 
5.'6 and 3.'2 to maximize the included emission. The spec- 
trum of NGC 6215A was extracted from a circular region 
of radius 4.'3. Source-free regions of each detector were 



used for the corresponding background measurements. 

With the limited data, only simple models are appropri- 
ate. In both cases, we fit the data reasonably well with a 
single power law. The model parameters are listed in Table 
3. NGC 6215 has been identified as a weak LINER (Durret 
& Bergeron 1988). The relatively strong X-ray emission of 
NGC 6215A suggests that this galaxy also contains an ac- 
tive nucleus. In addition to NGC 6215 and NGC 6215A, 
we observe the Seyfert 2 ESQ 138-Gl in the GIS fields. 
CoUinger & Brandt (2000) discuss the spectrum of ESQ 
138-Gl and find that either a partial covering or Compton 
reflection model fits it well. 

The GIS3 calibration source blocks the cluster 
1RXSJ165259.4-594302, and we observe it in the GIS2 
detector alone. The cluster is centered outside the GIS 
field, and we extract its spectrum from a 6^7 circular re- 
gion on the detector. We fit the spectrum with a thermal 
bremsstrahlung model having kT = 2.8 keV absorbed by 
fixed Galactic column density. The model parameters are 
listed in Table 3. 

5. THE OBSCURING STARBURST MODEL 

To account for the optical and X-ray observations of 
NGC 6221, we suggest a scenario in which the gas and dust 
associated with the starburst is the material that obscures 
the AGN. The extent of the star burst-like optical emission 
line ratios indicates that the physical scale of the starburst 
is a few hundred parsccs. This size is sufficient to block not 
only the broad-line region, but also the narrow-line region 
of the AGN. This is unlike the putative obscuring "torus" 
that distinguishes type 1 and type 2 Seyfert galaxies be- 
cause the latter hides only the broad-line region (Figure 
6). The total obscuring column density to the AGN mea- 
sured at X-ray energies is only 10^^ cm"^. The Seyfert 
1-like characteristics of NGC 6221, such as variability and 
Fe Ka line properties, are therefore not surprising. With- 
out the obscuring starburst, we would have a direct view of 
the AGN (Figure 6, viewpoint A). The column density to 
the nucleus of a true Seyfert 2 with an obscuring starburst 
is even greater (viewpoint B), as both the star- forming re- 
gion and the material very close to the central engine (the 
"torus") block the view (Levenson et al. 2001a). 
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Fig. 6. Cartoon of the obscuring starburst model. Wc view the nucleus of NGC 6221 from position A, through the gas and dust of the 
starburst. Without the starburst, NGC 6221 would appear to be a normal Seyfert 1. It is not a true Seyfert 2/starburst composite, in which 
our line of sight would pass through both the star-forming region and the "torus" very close to the central engine (viewpoint B). 



The starburst itself is dusty, as the optical spectra 
demonstrate. The extinction Ay = 3.0 mag (§2) is con- 
sistent with the model geometry, which places the star- 
burst along the line of sight to the AGN. For Nh = 
2 X lO^^^ycm-2 (Gorenstein 1975), we have Nh = 6 x 
lO^-'^cm"^. This measures the minimum obscuration due 
to the starburst, since most of the detected emission comes 
from its exterior regions that are not so extremely ab- 
sorbed. Specifically, the direct conversion from Ay to Nh 
assumes that all the dust is located in a uniform foreground 
screen of material, separate from the emission sources. As 
more complex model calculations demonstrate, the effec- 
tive extinction that we measure can severely underesti- 
mate the total dust and gas content of the starburst. In 
the starburst model of Witt et al. (1992), for example, in 
which the stars are more centrally c;oncentratcd than the 
embedding dust, the effective extinction is only Ay = 1.7 
even when the total extinction to the center of the galaxy 
reaches Av,totai = 15. Similarly, with a clumpy obscur- 
ing medium near the stars, Calzetti, Kinney, & Storchi- 
Bergmann (1994) find Ay = 2 would be measmed from 
the Balmer decrement when the total extinction along the 
line of sight Av,totai = 10, assuming a large number of 
individual clumps. 

In addition to the distinct physical location of the X-ray 
and optical emission, a second uncertainty in comparing 
the X-ray and optical measurements arises because they 



are directly sensitive to different obscuring material. Pho- 
toelectric interaction with heavy elements is the immediate 
culprit at X-ray energies, and determining the hydrogen 
column density then is a function of abundance. Simi- 
larly, the optical extinction is due to dust, and the equiva- 
lent hydrogen column density depends on the gas-to-dust 
ratio. Nevertheless, the lower limit on column density of 
the starburst is at least comparable to Nh that obscures 
the AGN as measured in X-rays, and more complex ge- 
ometries of emission sources and dust within the starburst 
indicate that it is likely to be greater, which implies that 
the starburst intervening along the line of sight to the ac- 
tive nucleus of NGC 6221 is the primary source of obscu- 
ration. Note that because the extinction we measure is to 
the extended starburst region, we expect standard Galac- 
tic Ay /Nh ratios to be appropriate, in contrast with some 
AGNs, in which extinction to the broad-line region is lower 
than expected for the absorbing column density measured 
in X-rays and may be due to fundamentally different dis- 
tributions of grain size or dust composition close to the 
central engine (e.g., Maiolino et al. 2001). 

The presence of both a starburst and an AGN in this 
geometry naturally results in observed [O III] lines that 
are broader than other lines if the two components are 
not resolved. The detectable AGN contribution to Balmer 
lines and other features that are characteristically strong 
in starbursts is relatively small. The starburst is relatively 
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weaker in [O III], however, so the fraction of observed lu- 
minosity in this hne that is due to the AGN is greater. The 
AGN Une broadening then weights the total measured line 
width more heavily. Broadened [O III] lines in fact in- 
dicate the presence of AGNs in substantial fractions of 
galaxies with starburst and H Il-region-like nuclei (Kenni- 
cutt, Keel, & Blaha 1989), in addition to NGC 6221 and 
the other X-ray loud composites. 

High spatial resolution observations of starburst galaxies 
demonstrate that the starbursts are inhomogenous (e.g., 
Calzetti et al. 1997), so in this case, we do not expect the 
AGN to be covered uniformly. While the physical scale of 
individual clumps may be as small as distinct star-forming 
clouds, however, the net obscuration to the center of NGC 
6221 is large. Similar to other starburst galaxies, in which 
total nuclear column densities of 10^^ to lO^^cm"^ are 
typical (e.g., Dahlem, Weaver, & Heckman 1998; Kenni- 
cutt 1998b), the line of sight intersects many dusty re- 
gions, so the AGN is never directly exposed. If changes 
in the intrinsic AGN and its obscuring medium can be 
distinguished. X-ray variability measurements may best 
constrain the physical scale of the obscuring clouds. 
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Fig. 7. — Net observed ratio, R, of [O III]/H/3 intensity vs. 
AGN fraction of total intrinsic H/3 flux, calculated for Ragn = 10, 
j^AGN _ g 2^ 4^ and 6 (increasing from upper left to lower right), 
and RsB = 0.1 and 1.0. With some obscuration, even a strong 
AGN may not be detected in the observed line ratios. The AGN 
fraction of observed H/S flux is marked with ticks corresponding to 
0.1, 0.2, 0.3, 0.4, and 0.5 along the Rsb = 0.1 curves. The rela- 
tionship between the observed and intrinsic H/3 fractions depends 
only on A^'^^, not Rsb or Ragn- For example, with A^'^^ = 2, 
the observed H/3 fraction equals 0.1 when the intrinsic Hp frax;tion 
is 0.4. 

The HST images demonstrate that the observed opti- 
cal flux due to the AGN is small compared with the sur- 
rounding starburst. Both sources are strongly extincted, 
and the AGN is also intrinsically weak with respect to 
the starburst at optical wavelengths. We roughly estimate 
the AGN properties from the X-ray luminosity. First we 
estimate the broad H/3 line luminosity from the hard X- 
ray luminosity, L2~io- Xu, Livio, & Baum (1999) compile 
the X-ray measurements, and we find the corresponding 
H/3 values of the Seyfert Is in Whittle (1992) and Da- 
hari & de Robertis (1988). The median ratio of these is 
Lh0^/L2-io = 0.013. In the high-flux state of NGC 6221, 



we predict if^^ = 8.5 x lO^'^ergs'S or flux = 

1.8 X 10^^'^ ergcm^^ s^-'^ from the intrinsic AGN. Expect- 
ing Ay = 5.0 to the broad-line region based on the X- 
ray column density, the resulting expected detectable flux 
F§^^ = 1.8 X 10"^^ ergcm-2 s-\ which is well below the 
upper limit on detection of the BLR estimated from the 
optical spectrum. Even adopting the minimum Ay =3.0 
to the AGN based on the optical data, the expected mea- 
sured flux F^^^ = 1.1 X 10~"ergcm-2s~\ at the detec- 
tion limit in the optical spectrum and less than the total 
observed Fnp ~ 2.4x10^^^ erg cm~^ s~^ in the nuclear re- 
gion. Using the X-ray-Balmer line luminosity correlation 
of Ho et al. (2001) for broad-line galaxies produces simi- 
lar results, with the weak BLR undectable in our optical 
spectrum. 
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Fig. 8. — Net observed ratio, R, of [O III]/H^ intensity vs. 
AGN fra<;tion of total intrinsic H/3 flux, calculated for Rsb = 0.1, 
^AGN _ g^ 2^ 4^ and 6 (increasing from upper left to lower right), 
and Ragn = 3 and 10. The AGN fraction of observed H/3 flux 
is marked with ticks corresponding to 0.1, 0.2, 0.3, 0.4, and 0.5 
along the Ragn = 3 curves. The combination of obscuration and 
relatively low AGN power can produce line ratios typical of pure 
starbursts. 

Second, we crudely estimate the intrinsic AGN 
[O III]A5007 luminosity, which is produced in the NLR, 
from ^2-10- The median ratio of these two quantities 
^[o/?/]/^2-io = 8.7 X 10-3 for the Seyfert Is of the 
Xu et al. (1999) compilation. For the observed L2-10 = 



6.6 X 10 erg cm 
or intrinsic flux F^^^^ 



5.7 X lO'^^ergcm ^s 



1 tNLR 
' ^[OIII] 

-- 1.2 X 10-"ergcm-2s-^. The 
advantage of calibrating this relationship with only the 

Seyfert Is is that the ratio L^^jfj^/ L2-1Q is sensitive 
to intervening column density, although fitting all ra- 
dio quiet galaxies in the data set produces similar re- 
sults. From the AGN component observed at optical 
wavelengths, wc expect the reddening-corrected J'jof/}] = 

4.2 X lO-^'^ergcm^^ s^^ in the 1'.'5 x 1'.'5 central aperture, 
or F^jpj-^ = 8.1 X lO^^'' erg cm"^ s"^, integrated over the 
whole slit. Two uncertainties relevant to both of these X- 
ray estimates are that the data set of X-ray and optical 
fluxes is not uniform and the optical fluxes of this sample 
have not been corrected for reddening. The scatter in the 
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data is large, which may therefore account for our undcr- 
prediction of the observed [O III] luminosity. Although the 
flux measured in the optical spectrum is approximately 
four times greater than the intrinsic [O III] flux we pre- 
dict based on the hard X-ray flux, our main point is to 
demonstrate that the X-ray and optical data all provide a 
generally consistent description of the starburst and AGN 
in NGC 6221. 
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Fig. 9. — Spectral energy distributions of unobscured sources 
illustrate the X-ray prominence of AGNs. The average intrinsic 
starburst {solid line; Schmitt et al. 1997) and radio-quiet quasar 
{dotted line; Elvis et al. 1994), are each normalized at 5000A. We 
model the X-ray spectrum of the starburst as thermal emission due 
to gas at kT = 1 keV, and the dashed line covers the UV region 
where there intrinsic spectrum is unknown. When the optical fluxes 
are equal, the AGN is two orders of magnitude brighter than the 
starburst in X-rays, although the bolometric luminosities of the dif- 
ferent emission sources are comparable. 

Mixing starburst and AGN components can produce an 

ambiguous optical spectrum, as Hill ct al. (2001) and Nel- 
son & Whittle (2001) demonstrate for several diagnostic 
line ratios. Additional obscuration of the AGN can com- 
pletely hide this component to produce an optical spec- 
trum that mimics a pure starburst. We compute some 
simple models to illustrate the detected changes in the 
important flux ratio R = F^oiii]/ Fnp in general circum- 
stances. In the first example (Figure 7), we adopt a typical 
AGN ratio Ragn = 10 and optical extinction to the star- 
burst Ay^ ~ 0, and we consider the effects of additional 
extinction to the AGN {A^'^^ = 0, 2, 4, and 6) and dif- 
ferent starburst line ratios {Rsb =0.1, and 1.0). Even 
when the intrinsic AGN is strong, providing up to 50% of 
the completely unobscured narrow H/3 luminosity, the net 
observed R can be indistinguishable from a pure starburst. 
Similarly, we consider Ragn = 3 and 10 and show that 
for the typical Rsb =0.1, some combination of obscura- 
tion and relative AGN luminosity can hide the AGN in the 
observed line ratio (Figure 8). In all cases, only the addi- 
tional obscuration of the AGN, Ay'^'^ — Ay^ , is relevant, 
and the plots are identical for A^^ = 3 and Ay'^^ = 3, 
5, 7, and 9, which encompasses the likely parameters in 
the case of NGC 6221. We also note the observed frac- 
tion of II/3 flux due to the AGN with ticks corresonding 
to observed fractions 0.1, 0.2, 0.3, 0.4, and 0.5. For ex- 
ample, with Ay*^^ = 2, 10% of the detected H/3 emission 



is due to the AGN when its intrinsic fraction is 0.4, while 
for Ay'^^ = 4, 10% of the detected H/3 emission is due to 
the AGN only when the an intrinsic AGN fraction reaches 
0.8. The relationship between the intrinsic and observed 
H/3 fractions depends only on Ay'^^ , not Rsb or Ragn- 
The effect of obscuration on typical spectral energy dis- 
tributions also illustrates that an AGN can be hidden at 
optical wavelengths while remaining the dominant X-ray 
source in composite galaxies. We compare the average un- 
obscured starburst galaxy (Schmitt ct al. 1997) with the 
median radio-quiet QSO (Elvis et al. 1994), normalized at 
5000A (Figure 9). The bolometric luminosities of the two 
sources are comparable, with the starburst much stronger 
in the FIR and the AGN dominating in X-rays. 



Icm Imm 100|J.m lO^m l|im lOOOA O.lkeV IkeV 



44 


1 1 1 1 


=5x10" 


cm' 






N„"»'=lxlO=cm-' 


42 












> 










40 








38 








36 









10 12 14 16 18 

log V (Hz) 

Fig. 10. — Observed spectral energy distribution of NGC 6221, 
modelled with typical AGN and starburst contributions. The AGN 
component is measured in X-rays and as the point source alone in 
B.ST images {filled symbols). Imaging and spectroscopy of NGC 
6221 that encompass larger scales are predominantly sensitive to 
the starburst (open symbols). The average AGN (dotted line) is ex- 
tincted by 1 X 10^^ cm~^, as we measure in the X-ray spectrum of 
NGC 6221, while the average starburst (solid line) is cxtincted by 
5 X 10^1 cm-^, as we measure optically. In NGC 6221, the AGN is 
intrinsically weak compared with the starburst, as well as heavily 
absorbed. While the starburst accounts for nearly all the bolo- 
metric luminosity of NGC 6221, the AGN still produces its X-ray 
emission. Additional data in the FIR, optical-NIR, and UV bands 
are from Sanders et al. (1995), McAlary et al. (1983), and Kinney 
et al. (1993), respectively. 

We model NGC 6221 with these two components, scal- 
ing the empirical spectral energy distributions to the data 
values and applying the measured extinctions (Figure 10). 
The X-ray and HST measurements of the point source in- 
dicate the contribution of the AGN, and we use the X-ray 
determination of its extinction. Lower-resolution imag- 
ing and large-aperture spectroscopy measure the starburst 
component, and the Balmer decrement indicates its exct- 
inction. In this calculation, we adopt the extinction laws of 
Cardelli ct al. (1989) at long wavelengths, that of Calzctti 
et al. (2000) in the NlR-to-UV bandpass, and that of Mor- 
rison & McCammon (1983) at higher energies. While the 
empirical Calzetti et al. (2000) law is based on observa- 
tions of starburst galaxies, it is also reasonably applied to 
the AGN component. Specifically, it is relatively grey and 
lacks the strong Galactic feature at 2 175 A, which is not 
detected in distant quasars (Pitman, Clayton, & Gordon 
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2000, and references therein). The AGN is much weaker 
than the starburst at aU wavelengths except X-rays. The 
AGN accounts for only about 1% of the galaxy's bolo- 
metric luminosity, while providing nearly all the observed 
X-rays. In NGC 6221 and other X-ray loud composites, 
the dusty starburst is both the obscuring medium and the 
origin of the detectable optical and infrared spectrum. 

The population of starburst-obscured AGNs helps to 
solve the problem of the spectrum of the XRB. The ob- 
served XRB spectrum is miich flatter than the spectra of 
unobscured AGNs or any other extragalactic populations 
that likely comprise the background. A combination of ab- 
sorption and rcdshift of AGNs can produce the observed 
XRB, however (Setti & Woltjer 1989; Madau, Ghisellini, 
& Fabian 1994; Comastri ct al. 1995). As Fabian et al. 
(1998) have proposed, starburst-obscured AGNs specifi- 
cally may contribute significantly to the XRB, naturally 
providing the required large optical depths and high cov- 
ering fractions of their central engines. 

The starburst and AGN phenomena may be related in 
an evolutionary sequence, with galactic interactions trig- 
gering the starbursts and providing a central reservoir of 
material for accretion in the active nucleus (Cid Fcrnan- 
des et al. 2001; Storchi-Bergmann et al. 2001). The lim- 
ited lifetime of the starburst then constrains the timescale 
of the strongly obscured AGN phase. Typical starbursts 
persist only over timescales up to several hundred Myr. 
As a result of the starburst, material is blown out of the 
galaxy by the superwind or converted into dense stellar 
remnants, and no longer blocks the active nucleus. Thus, 
as Heckman et al. (1989) describe, only in a later phase 
of evolution are direct views of the central engine possi- 
ble, with characteristic optical spectra that include broad 
emission lines. 

6. SUMMARY AND SIGNIFICANCE 

We have presented NGC 6221 as a case study of an X-ray 
loud composite galaxy. While the optical spectrum of the 
nuclear region is typical of a starburst galaxy, its buried 
AGN is clearly revealed in X-rays. The Balmer decrement 
indicates that the starburst itself is strongly reddened, and 
the column density to the AGN measured in X-rays is even 
greater. We suggest that the very red point source distin- 
guished in HST images is the nucleus, and we crudely de- 
rive intrinsic optical properties of the AGN from the X-ray 
luminosity. We demonstrate generally that the combina- 
tion of obscuration and contrast in viewing an intrinsically 
weak AGN against a bright starburst can produce the net 
optical spectrum of the surrounding starburst, while the 
active nucleus remains the dominant X-ray source. 

Without the obscuring starburst, NGC 6221 woidd be 
classified as a Seyfert 1, with the central engine and broad- 
line region viewed directly. In the simple physical model, 
the starburst need not completely block the AGN from all 
lines of sight, but its covering fraction may be larger than 
that of the dense material very close to the nucleus — the 



"torus" — that accounts for some of the observed distinc- 
tions between Seyfert Is and 2s. In some cases, the star- 
burst may block only the broad emission lines, not the 
narrow ones. Consequently, if such obscuring starbursts 
arc common, the observed ratios of Seyfert Is and 2s do 
not directly constrain the "torus" characteristics such as 
covering fraction. 

NGC 6221 is not unique, and such objects probably oc- 
cur more frequently than their small fraction (7/210) in 
the joint RO S AT- IR A S s-Axn^Xc (Boiler ct al. 1992; Moran 
et al. 1996) would suggest. Because the ROSAT survey 
is sensitive only to soft (0.2-2 keV) X-rays, it misses the 
more heavily obscured members of this class. A column 
density of 10^^ cm~^ reduces the flux in the ROSATha,nd- 
pass by a factor of 10, and Nh = 10^'^ cm~^ reduces the 
flux by more than three orders of magnitude. These col- 
umn densities are not extreme compared with the ob- 
served nuclear column densities of molecular gas in or- 
dinary starburst galaxies. While optical spectroscopy is 
the primary method of identifying the sources detected 
through efforts to resolve the X-ray background, it does 
not always reveal the complex nature of these galaxies. 
Inverting the problem - considering the X-ray emission of 
known optically-identified populations — also overlooks the 
contribution of the X-ray loud composite AGNs. Hard 
X-ray surveys are the best method to detect this class of 
galaxies, and they require careful spectroscopy at lower en- 
ergies to correctly identify their intrinsic; emission sources. 
Future missions with greater sensitivity to harder X-rays 
will prove most effective at finding X-ray loud AGNs. 
Constellation- X, for example, is expected to have an ef- 
fective area of about 5000 cm^ at 10 keV and 1500 cm^ 
at 40 keV, and will therefore sensitively detect obscured 
AGNs of all kinds. 
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